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Abstract 

This  work  has  two  parts.  The  first  part  details  the  thermal  analysis  of  a  thin-film  thermoelectric  cooler  under  the  influence  of  the  Thomson  heating, 
the  Joule  heating,  and  the  Fourier’s  heat  conduction.  A  constant  Thomson  coefficient,  instead  of  traditionally  a  constant  Seebeck  coefficient,  is 
assumed.  The  influence  of  the  Thomson  effect  on  the  cooling  power,  the  achievable  temperature  difference  and  the  optimum  operating  current 
density  is  then  explored.  It  is  found  that  the  Joule’s  heat  and  the  conduction  heat  flowing  to  the  cold  junction  can  be  significantly  reduced  if 
thermoelectric  materials  with  properly  designed  Thomson  coefficients  are  employed.  A  modified  thermal  conductance  and  a  modified  electric 
resistance  are  resulted. 

The  second  part  of  this  paper  details  the  analysis  of  the  thermal  stresses  existing  in  the  layered  structure  and  induced  by  the  temperature  difference 
via  a  non-coupled  thermal  elastic  theory.  The  results  provide  a  preliminary  knowledge  to  judge  whether  the  thin-film  structure  is  destroyed  by  the 
thermal  stresses  or  not,  especially  by  the  shear  stresses  between  adjacent  layers. 

©  2005  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Thermoelectric  cooler;  Thomson  effect;  Thermal  analysis;  Thermal  stress 


1.  Introduction 

Thermoelectric  coolers  have  advantages  over  the  traditional 
heat  pumps  because  they  are  compact  and  highly  reliable,  have 
no  moving  parts,  and  use  no  refrigerants.  The  poor  cooling 
efficiency,  which  is  fundamentally  limited  by  the  material  prop¬ 
erties,  of  thermoelectric  devices,  however,  limits  their  applica¬ 
tions.  Basically,  the  cooling  power  arises  from  the  Peltier  effect 
and  is  degraded  by  the  Fourier’s  conduction  heat  due  to  the 
existing  temperature  difference  and  by  the  Joule’s  heat  due  to 
the  electric  resistance.  A  material  with  a  large  Seebeck  coeffi¬ 
cient  (a),  a  large  electrical  conductivity  (er),  and  a  small  thermal 
conductivity  (k)  is  thus  preferred.  In  fact,  a  conventional  ther¬ 
mal  analysis,  taking  the  Fourier’s  conduction  heat  transfer,  the 
Joule’s  heating,  and  sometimes  the  radiation  and  convection  heat 
transfer  between  the  thermoelectric  element  and  the  ambient 
gas  into  consideration  [1-3],  shows  that  thermoelectric  materi¬ 
als  are  estimated  by  their  ZT  values  (figure-of-merit),  where  T  is 
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the  absolute  temperature  and  Z=aa2lk.  The  commonly  known 
thermoelectric  materials  have  ZT  values  between  0.6  and  1.0  at 
room  temperature.  It  is  believed  that  practical  applications  could 
be  many  more  if  materials  with  ZT  values  greater  than  3  could 
be  developed. 

Various  efforts  [4-7]  have  been  made  to  develop  materials 
with  higher  ZT  values.  Recently,  some  exciting  results  have  been 
reported  at  the  material  research  society  meetings  for  ZT>  2  in 
use  of  some  low-dimensional  materials  such  as  quantum  wells, 
quantum  wires,  quantum  dots,  and  superlattice  structures  [8-15] . 
The  increase  in  the  ZT  values  is  explained  by  the  belief  that 
reduced  dimensionality  changes  the  band  structures  (enhances 
the  density  of  states  near  the  Fermi  energy),  modifies  the  phonon 
dispersion  relation,  and  increases  the  interface  scattering  of 
phonons.  Consequently,  the  electric  resistance  and  the  lattice 
thermal  conductivity  [16-19]  are  both  reduced,  particularly  the 
latter.  In  labs,  thin-film/ superlattices  thermoelectric  devices  with 
very  small  dimensions  have  been  fabricated  using  microelec¬ 
tronics  technology  and  quantum  wires  are  in  fabricating. 

A  careful  examination  on  the  conventional  thermal  analysis 
finds  that  the  temperature  dependence  of  the  Seebeck  coefficient 
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is  usually  ignored  and  therefore  the  Thomson  effect  (charge 
carriers  must  absorb  or  release  heat  in  order  to  overcome  the 
non-homogeneity  caused  by  the  temperature  gradients  in  the 
semiconductor)  is  not  considered  [20,21],  Nonetheless,  the  tem¬ 
perature  dependence  of  the  Seebeck  coefficient  is  often  observed 
in  the  labs  [22-24]  and  the  temperature  gradients  exist  surely  in 
a  thermoelectric  device.  Hsu  et  al.  [22]  measured  the  thermo¬ 
electric  properties  of  bulk  AgPbioSbTei2-  From  their  Fig.  2A, 
the  Thomson  coefficient  is  estimated  to  be  about  the  same  as 
the  Seebeck  coefficient  («^200  p,V/K).  The  temperature  depen¬ 
dence  of  bismuth  films  of  different  thicknesses  was  explored 
by  Damodara  and  Soundararajan  [23].  Their  measurements  also 
showed  a  Thomson  coefficient  on  the  order  in  magnitude  of  the 
Seebeck  coefficient.  The  signs  are  the  same  near  the  room  tem¬ 
perature  but  opposite  at  higher  temperatures.  Experimental  data 
obtained  by  Rowe  et  al.  [24]  for  intermediate  valence  compound 
YbAl3  [24]  shows  similar  results. 

Not  impossible  is  that  the  reversible  Thomson  heat  changes 
the  temperature  profile  in  the  thermoelectric  element  and  modi¬ 
fies  the  associated  thermal  conductance  and  electric  resistance, 
resulting  in  a  better  thermal  performance  [25,26].  This  con¬ 
stitutes  the  motivation  of  the  present  study,  to  investigate  the 
influence  of  the  Thomson  effect  on  the  performance  of  a  thin- 
film  thermoelectric  cooler.  On  the  other  hand,  although  a  better 
thermal  performance  is  targeted,  the  thermal  stresses  induced  by 
the  raised  temperature  difference  may  be  harmful  [27].  Above 
all,  the  induced  shear  stresses  between  adjacent  thin-film  layers 
may  tear  the  layered  structure.  A  one-dimensional  thermal  stress 
analysis  therefore  will  also  be  performed. 

This  paper  is  arranged  as  follows.  The  thermal  model  adopted 
is  described  and  the  thermal  performance  in  the  presence  of  the 
Thomson  effect  is  investigated  in  Section  2.  The  thermal  stress 
analyses,  including  the  normal  stresses  and  the  shear  stresses, 
are  performed  through  a  non-couple  elastic  theory  and  presented 
in  Section  3.  Conclusions  are  given  at  last  in  Section  4. 

2.  Thermal  analysis 

2.7.  Thermal  model 


radiation  and  convection  heat  transfer  between  the  thermocouple 
and  the  ambient  gas  is  ignored  in  the  present  study.  All  the  ther¬ 
moelectric  properties  of  the  employed  materials  are  assumed  to 
be  constant,  except  the  Seebeck  coefficient.  The  Seebeck  coef¬ 
ficient  (a)  instead  is  assumed  to  have  a  logarithmic  dependence 
on  the  temperature,  that  is 

«  =  «0  =  £ln  J  (1) 

T0 

where  ao  is  the  Seebeck  coefficient  at  some  reference  tempera¬ 
ture  To.  In  other  words,  it  is  the  Thomson  coefficient,  ft,  that  is 
assumed  to  be  temperature-independent  [20]. 

According  to  the  non-equilibrium  thermodynamics  and  fol¬ 
lowing  the  analysis  of  Volklein  et  al.  [3],  one  can  write  down  the 
governing  equation  of  the  steady  temperature  distribution  T(x) 
in  the  thermocouple  as  follows: 


°-iK  Wau 


(2) 


where  I  is  the  operating  electric  current,  flowing  from  x  =  0  to 
x  =  L  in  the  p-type  semiconductor  and  from  x  =  L  to  x  =  0  in  the 
n-type  semiconductor;  /Jpn  =  f3p  —  and 


(kA\  =  kpAp  +  kiAiknAn  +  kmAm 

1  _  1  1 

(n/t)pn  ~  o-pAp  <rnAn 


(3) 

(4) 


The  subscripts  m,  p,  n,  and  i  represent  the  membrane,  the  p- 
type,  the  n-type,  and  the  insulating  layer,  respectively.  Finally, 
A  is  the  cross-sectional  area  of  each  layer.  Terms  on  the  right- 
hand-side  of  Eq.  (2)  are  the  net  axial  Fourier’s  conduction  heat 
transfer  rate,  the  Joule’s  heat  generation  rate,  and  the  Thomson 
heat  generation  rate,  respectively. 

With  the  boundary  conditions  7(0)  =  Tc  (the  cold-junction 
temperature)  and  T(L)  =  7),  (the  hot-junction  temperature),  the 
solution  of  Eq.  (2)  is  found  to  be 


T(x)  -  Tc 
Th~Tc 


=  (i~0 


1  ~  exp  (Hj) 

1  -  exp© 


(5) 


Considered  herein  is  a  thin-film  thermocouple  similar  to  that 
analyzed  by  Volklein  et  al.  in  1999  [3].  Fig.  1  shows  the  cross- 
section  of  the  structure,  including  the  supporting  membrane  and 
the  p-type  and  n-type  thermoelectric  films  that  are  separated  by 
an  insulating  layer.  As  mentioned  by  Volklein  et  al.  [3],  because 
the  thickness  of  the  membrane,  of  the  thermoelectric  films,  and 
of  the  insulating  layer  are  very  thin  compared  to  the  length  L,  the 
temperature  variations  in  the  thickness  direction  are  expected  to 
be  negligible.  Therefore,  a  one-dimensional  thermal  analysis  is 
adopted.  Moreover,  in  order  to  highlight  the  Thomson  effect,  the 


I  I  p-type 

insulator 
n-type 
membrane 

Fig.  1 .  The  cross-section  of  the  thin-film  thermocouple  under  investigation. 


where  f  =  /tpn//k't  and  £  =  //?t//3pn(Th— Tc)  are  the  ratio  of  the 
associated  Thomson  heat  to  the  Fourier’s  conduction  heat  and 
the  ratio  of  the  Joule’s  heat  to  the  Thomson  heat.  The  “intrinsic” 
thermal  conductance  and  electric  resistance  of  the  thermocouple 
(those  in  the  absence  of  the  Thomson  effect)  are  designated  as 
Kt  =  (kA)t/L  and  Rt  =  U(aA)pn.  When  the  Thomson  effect  is 
absent,  Eq.  (5)  is  reduced  to 


T(x)  -  Tc 
Th~Tc 


(6) 


in  which  the  ratio  of  the  Joule’s  heat  to  the  Fourier’s  conduction 
heat  (=££)  alone  determines  the  temperature  distribution.  The 
influence  of  the  Thomson  heat  on  the  temperature  profiles  is 
shown  in  Fig.  2,  in  which  ££  is  fixed  to  be  2  such  that  the  dimen¬ 
sionless  temperature  distribution,  Eq.  (6),  has  a  maximum  value 
of  1  at  x  =  L.  An  increasing  temperature  gradient  near  x  =  L  and 
a  decreasing  one  near  x=0  with  an  increasing  positive  f  are 
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x/L 


Fig.  2.  The  dimensionless  temperature  profiles  within  the  thermoelectric  ele¬ 
ment  under  several  ratios  of  the  Thomson  heat  to  the  Fourier’s  conduction  heat 
(f)  with  ff  =  2  fixed. 


under  a  given  temperature  difference  AT  =  7h  —  Tc  and  an  oper¬ 
ating  current  I  is  computed  as 

N  =  apnITc  -  K(Th  -  Te)  -  ~RI2/2  (10) 

2.2.  Thermal  performance 

The  performance  of  a  thermoelectric  cooler  is  judged  by  its 
maximum  achievable  temperature  difference  A  Tmax  or  by  its 
maximum  allowable  heat  load  Nmax.  The  former  is  the  maximum 
temperature  difference  that  can  be  achieved  under  no  external 
heat  load  as  the  operating  current  varies.  From  Eq.  (10),  it  can 
be  shown  that  the  maximum  achievable  temperature  difference 
A Tmax  and  the  corresponding  optimum  operating  current  Iopt.T- 
are  determined  by 

ZCTC  =  ^(1  -  e“f°PtT)  (11) 

Ppn 

and 


observed.  The  latter  implies  an  improved  thermal  performance 
of  the  thermoelectric  cooler  by  the  Thomson  effect. 

Based  on  the  temperature  distribution,  Eq.  (5),  the  Fourier’s 
heat  flow  rate  at  the  cold  junction  can  be  calculated  and  written 
as 

(kA\  ~  =  K(Th  -  Tc)  +  RI2/ 2  (7) 

dx  L=o 

where  K  =  and  R  =  rjR(^)Rt  are  the  modified  (effec¬ 

tive)  thermal  conductance  and  the  modified  (effective)  electric 
resistance.  The  modification  factors  are 

(8) 

vr(S>  =  i  -  -F— r  (9) 

§  e«  —  1 

and  are  shown  in  Fig.  3.  It  is  interesting  to  see  that  both  decrease 
with  increasing  positive  §pn.  Finally,  the  allowable  heat  load 


“nn 

ZcATmax  =  ^  {f0pt.T  +  exp(— £optT)  -  1}  (12) 

where  £0ptT  =  PpJopijIKu  apn  =  ap  -  an  and  Zc  =  a^JKJRt 
(evaluated  at  the  cold-junction  temperature  Tc).  For  a  compar¬ 
ison,  the  counterparts  when  the  Thomson  effect  is  absent  are 
written  below: 


/ opt.T  — 
ATjnax  = 


OjmTc 


(13) 

(14) 


Fig.  4  shows  one  illustration.  The  material  properties  employed 
by  Volkleinetal.  [3]  are  referred.  They  are:  rrf]n  =  13.3  |x£2  /m, 
kpn  =  1 .28 W/mK,  lijm  =  2.0W/mK,  and  o:pn  =  370 p,V/K  at 
Tc  =  250  K.  The  thermocouple  dimensions  are  chosen  to 
be  Aj+m  =  0.2  mm  x  3  pm  x  2,  Ap  n  =  0.2  mm  x  10  pm,  and 
L  =  0.3  mm.  As  seen,  the  maximum  achievable  temperature  dif¬ 
ference  is  increased  by  a  positive  fpn.  It  is  about  42.8  K  when 
Ppa  =  0  and  about  57.5  K  when  /3pn  =  400  pV/K.  This  is  because 


I 


P(gV/JO 


Fig.  3.  The  modification  factors  of  the  thermal  conductance  and  the  electric  Fig-  4-  The  influence  of  Thomson  effect  on  the  maximum  achievable  tempera- 
resistance  due  to  the  Thomson  effect.  ture  difference  and  the  corresponding  optimum  operating  current  density. 
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the  current  absorbs  heat  not  only  at  the  cold-junction  through 
the  Peltier  effect  but  also  all  the  way  along  the  thermoelectric 
element  due  to  the  Thomson  effect.  The  amounts  of  Joule’s  heat 
and  Fourier’s  conduction  heat  flowing  to  the  cold  junction  are 
consequently  decreased  by  factors  of  or  and  ?7k,  respectively. 

The  maximum  allowable  heat  load,  /Vmax,  is  defined  as  the 
maximum  heat  load  under  AT  =0  as  the  operating  current  I 
varies.  Substituting  Eqs.  (8)  and  (9)  into  Eq.  (10),  one  finds 
that 

ZCTC  =  ^  ( 1  +  ^°pt  N  6Xp  ^opt'N  +  ZfoptN  \  ^ 

fipn  Y  (1  —  exp^0pt.N)2  1  —  exP^opt.N  J 

and 

^  yy  _  ^Pn  A  §opt.N  CXp  ^opt.N  \ 

t  PL  (1  -  exp^opt.N)  v  1  -  exp^opt.N  J 

(16) 


where  f  0pt.N  =  /fpn/opt.N/^t  is  the  corresponding  optimum  dimen¬ 
sionless  operating  current.  The  counterparts  as  /3pn  =  0  are 


(17) 


and 

AW  =  WtZcT2  (18) 

Note  that  the  two  optimum  operating  currents,  computed  from 
Eqs.  (11)  and  (15),  are  in  general  different.  Eqs.  (13)  and  (17) 
however,  are  the  same.  Fig.  5  shows  another  illustration.  The 
thermocouple  used  in  Fig.  4  is  employed  again  but  the  thermo¬ 
electric  properties  are  evaluated  at  Tc  =  7h  =  293  K  now  and  they 
are:  apn  =  460  p, V/K,  ct"*  =  17|xfi/m,  /fcpn  =  1 08  W/mK,  and 
AijjtT,  =  2.2  W/Mk.  The  maximum  allowable  heat  load  is  found  to 
be  17.8 mW  as  /3pn  =  0  and  is  19.4mW  when  /3pn  =  400  pV/K. 
One  may  notice  that  the  maximum  allowable  heat  load  is  not 
improved  as  much  as  the  maximum  achievable  temperature 
difference  is.  This  is  because  the  temperature  gradient  is  now 
positive  near  x  =  0  and  negative  near  x  =  L  (the  peak  tempera- 


and  the  corresponding  optimum  operating  current  density. 


x/L 


Fig.  6.  The  influence  of  Thomson  effect  on  the  temperature  profiles  when  the 
temperature  difference  across  the  thermoelectric  element  is  zero. 


ture  occurs  somewhere  in  the  middle  of  the  thermocouple).  The 
Thomson  effect  works  on  shifting  the  peak  temperature  toward 
the  hot  junction  and  reducing  the  peak  value  as  seen  from  Fig.  6, 
in  which  the  temperature  profiles,  Eqs.  (5)  and  (6),  are  re-plotted 
with  AT=  0.  As  a  result,  the  current  absorbs  heat  somewhere  but 
releases  heat  somewhere  else  along  the  thermocouple.  The  net 
improvement  is  thus  not  much. 

3.  Thermal  stress  analysis 

3.1.  Thermal  stress  model 

Like  the  thermal  analysis  performed  in  the  previous  section, 
it  is  also  assumed  that  the  thermal  displacement  occurs  mainly 
in  the  x-direction  and  its  variation  across  the  thickness  direction 
is  much  smaller  than  that  in  the  x-direction;  that  is,  u^>w  and 
9/9xi»  9 z.  The  free  stress  temperature  (7>s)  is  set  at  293  K  and 
non-elastic  effects  such  as  plasticity  and  creep  are  not  taken  into 
consideration.  The  thermoelectric  element  is  next  assumed  to 
be  clamped  at  both  junctions  and  free  at  the  top  and  bottom 
surfaces.  Consequently,  a  force  balance  over  a  control  volume 
as  shown  in  Fig.  7  results  in 

d^pOxx.p  d A\axx  i  dAncrxx  n  d Amcrxx  m  _ 

dx  +  dx  +  Ax  +  dx  “  1  ’ 

where  axx  is  the  normal  stress  and  is  related  to  the  x-displacement 
u(x)  by 

axx  =  E{jx~Y{T~ Tfs))  (20) 


Fig.  7.  The  control  volume  for  the  thermal  stress  analysis. 


126 


M.-J.  Huang  etal. /Sensors  and  Actuators  A  126  (2006)  122-128 


with  E  being  the  Young’s  modulus  and  y  being  the  thermal 
expansion  coefficient.  Substituting  Eq.  (20)  into  Eq.  (19),  one 
obtains 


d  2u  d  T 


(21) 


where 

_  ^p£pyp  +  AiEiVi  +  AnEnVn  +  AmEmym 

VO  A  1/tr^i/tTT  ,  . .  ("^i) 

ApEp  -4-  AiEi  -(-  AnEfi  -(-  AjyiEjyi 

is  the  weighted  average  of  the  thermal  expansion  coefficients. 
The  solution  of  Eq.  (21)  with  clamped  junctions  is  easily  found 
to  be 


«(*)  =  Kx_  /x  _  \  1  -K  /1-expgf 

Ly0AT  2  L\L  )  f  Vl-exp§ 


By  substituting  Eq.  (23)  into  Eq.  (20),  one  obtains  the  normal- 
stress  distribution  as  follows: 


Fig.  8.  The  dimensionless  displacement  along  the  thermoelectric  element  under 
several  ratios  of  the  Thomson  heat  to  the  Fourier’s  conduction  heat  (f)  with 
f  f  =  2  fixed. 


0«  =  E{yo(T  -  To)  -  y(T  -  Tfs)}  (24) 

where  the  current-dependent  temperature  7o  is  a  linear  combi¬ 
nation  of  the  junction  temperatures, 

T0  =  (pTh  +  (1  -  0)TC  (25) 


those  of  SiCE  are  used  for  the  membrane  and  the  insulating  thin- 
film  (£=160  Gpa  and  j/  =  5x  10-7/K).  As  seen,  with  increasing 
Thomson  effect,  the  absolute  normal  stresses  increase.  All  the 
calculated  normal  stresses  are  under  the  yielding  stresses  of  the 
materials  in  use  nonetheless. 


The  fraction  <p  is 


3.2.  Shear  stresses 


1  -  exp£  £ 
when  /J  ^  0  and 


1 


(26) 


(27) 


when  fi  =  0.  An  examination  on  Eqs.  (24)  and  (20)  shows  that 
the  temperature  To  appears  at  the  location  where  the  maximum 
displacement  occurs.  The  temperature  7o  may  be  viewed  as  a 
“free-structure-stress”  temperature.  Actually,  Eq.  (24)  implies 
that  the  thermal  stress  within  each  layer  can  be  divided  into  two 
parts.  The  first  part  is  caused  purely  by  the  local  temperature 
change  as  if  there  were  no  adjacent  layers.  The  second  part  on 
the  other  hand  is  induced  by  the  mutual  influence  between  layers 
due  to  the  layered  structure .  It  vanishes  if  the  temperature  is  equal 
to  To- 

Fig.  8  shows  the  influence  of  the  Thomson  effect  on  the  ther¬ 
mal  displacement.  It  is  observed  that  the  location  where  the 
maximum  displacement  or  7q  appears  shifts  toward  the  hot  junc¬ 
tion  as  the  Thomson  effect  increases.  The  dependence  of  the 
fraction  0  on  the  parameter  £  with  several  fixed  values  of  ££  is 
shown  in  Fig.  9.  Basically,  the  fraction  0,  or  (To—Tc)/(Th—Tc), 
decreases  with  increasing  Thomson  effect.  The  above  two  obser¬ 
vations  together  imply  a  reduced  temperature  gradient  at  the  cold 
junction  and  consequently  an  improved  thermal  performance. 

At  last,  presented  in  Fig.  10  are  the  normal-stress  distribu¬ 
tions  along  the  thermoelectric  element  operated  at  the  optimum 
current  shown  in  Fig.  5.  The  thermo-mechanical  properties  of 
semiconductor  Sio.sGeo.s  are  employed  for  the  p-type  as  well 
as  the  n-type  thin-film  (E=  160  Gpa  and  y  =  4.2  x  10-6/K)  and 


The  shear  stresses  between  adjacent  layers  can  be  found  now 
by  taking  a  force  balance  for  each  layer.  The  result  is 

dhcfrr 

rzx(x,  z  =  h)  =  —  +  zzx(x,  z  =  0)  (28) 

dx 

where  h  is  the  thickness  of  the  layer  of  interest.  Because  the 
top  and  the  bottom  surfaces  of  the  thermoelectric  element  are 
assumed  to  be  free,  the  shear  stresses  are  found  to  be 

dT 

Tm,n  =  -Emhm(yo  -  Ym)—  (29) 


4 


Fig.  9.  The  function  0  against  the  Thomson-effect  Parameter  f  when  the  value 
of  is  fixed  (the  increment  is  Af f  =  0.25). 
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(a)  x/L 


(b)  x/L 

Fig.  10.  The  normal-stress  distributions  of  the  thermoelectric  element  oper¬ 
ated  at  the  optimum  currents  shown  in  Fig.  5,  where  A/jpn  =  50  p,V/K  and 
y3pn,max  =400  pV/K.  The  dotted  curve  is  the  one  corresponding  to  /8pn  =  0. 


x/L 

Fig.  1 1 .  The  shear  stress  distributions  between  the  insulating  layer  and  the  p-type 
semiconductor  as  the  thermoelectric  cooler  is  operated  at  the  optimum  currents 
shown  in  Fig.  5,  where  A/3pll  =  50pV/K  and  ^pn,  max  =  400  pV/K.  The  dotted 
curve  is  the  one  corresponding  to  j6pn  =  0. 

distribution  between  the  p-type  layer  and  the  insulating  layer  is 
the  same  as  that  between  the  membrane  layer  and  the  n-type 
semiconductor.  Furthermore,  the  shear  stress  between  the  n- 
type  semiconductor  and  the  insulating  layer  vanishes.  Therefore, 
only  one  distribution  for  each  Thomson  coefficient  is  plotted  in 
Fig.  1 1 .  These  calculations  show  that  the  magnitudes  of  the  shear 
stresses  near  the  hot  junction  increase  more  rapidly  than  those 
near  the  cold  junction.  It  is  again  because  the  current  absorbs 
heat  near  the  cold  junction  and  releases  heat  near  the  hot  junction 
by  the  Thomson  effect,  causing  larger  temperature  gradients  (as 
observed  in  Fig.  6)  and  thus  larger  shear  stresses  near  the  hot 
junction. 

4.  Conclusion 


Tn.i  =  -[£nfcn(n>  -  Yn)  +  Emhm{yo  -  ym)]  ^  (30) 

and 

Ti.p  =  -[EMyo  -  Yi)  +  Enhn(yo  -  Yn) 
d  T 

+  Emhm(y0  -  Ym)l  .  (31) 

ax 

Alternatively,  the  shear  stress  between  the  insulating  layer  and 
the  p-type  layer  can  be  written  as 
dT 

Ti,p  =  Ephpiyo  -  yp)—  (32) 

due  to  Eq.  (19).  Eqs.  (29)-(32)  indicate  that  the  magnitude  of 
the  shear  stress  is  proportional  to  the  temperature  gradient.  In 
addition,  the  larger  the  heterogeneity  in  the  thermal  expansion 
coefficients  of  materials,  the  larger  the  shear  stress  is  resulted. 
From  Fig.  2,  the  location  where  the  maximum  temperature  gra¬ 
dient  appears  moves  gradually  from  the  cold  junction  to  the  hot 
junction  as  the  Thomson  effect  increases.  Fig.  1 1  shows  the  shear 
stress  distributions  corresponding  to  the  normal-stress  distribu¬ 
tions  in  Fig.  10.  Because  of  the  materials  used,  the  shear  stress 


The  thermal  performance  of  a  thin-film  thermoelectric  cooler 
under  the  influence  of  the  Thomson  heating,  the  Joule  heat¬ 
ing,  and  the  Fourier’s  heat  conduction  is  investigated.  For  those 
cases  that  the  Thomson  effect  is  important,  two  dimensionless 
influence  parameters — the  ratio  of  the  Thomson  heat  to  the  con¬ 
duction  heat  and  the  ratio  of  Joule’s  heat  to  the  Thomson  heat, 
determine  the  thermal  performance.  The  conduction  heat  and 
the  Joule’s  heat  flowing  to  the  cold  junction  can  be  significantly 
reduced,  providing  that  the  Thomson  coefficient  of  the  p-type 
layer  is  greater  than  that  of  the  n-type.  A  larger  maximum  tem¬ 
perature  difference  can  be  achieved  and  a  larger  maximum  heat 
load  can  be  allowed  consequently. 

The  normal-stress  distributions  within  the  thermoelectric  ele¬ 
ment  as  well  as  the  shear  stress  distributions  between  adjacent 
layers  are  also  analyzed  by  a  non-coupled  elastic  theory.  A  “ffee- 
structure-stress”  temperature  for  the  whole  thin-film  thermo¬ 
couple  is  defined.  It  characterizes  the  mutual  influence  between 
layers.  As  the  thermal  performance  is  improved  by  the  Thomson 
effect,  the  induced  thermal  stresses  are  increased  as  well.  They 
cannot  be  greater  than  the  yielding  stresses  of  the  materials  in  use 
for  safety.  Finally,  it  is  found  that  the  shear  stress  is  proportional 
to  the  temperature  gradient  and  the  larger  the  heterogeneity  in 
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the  thermal  expansion  coefficients  of  materials,  the  larger  the 
shear  stress  is  resulted. 
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